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Abstract 



We study one-magnon excitations in a random ferromagnetic Heisenberg 
chain with long-range correlations in the coupling constant distribution. By 
employing an exact diagonalization procedure, we compute the localization 
length of all one-magnon states within the band of allowed energies E. The 
random distribution of coupling constants was assumed to have a power spec- 
trum decaying as S(k) oc l/k a . We found that for a < 1, one-magnon excita- 
tions remain exponentially localized with the localization length £ diverging 
as 1/E. For a = 1 a faster divergence of £ is obtained. For any a > 1, a 
phase of delocalized magnons emerges at the bottom of the band. We charac- 
terize the scaling behavior of the localization length on all regimes and relate 
it with the scaling properties of the long-range correlated exchange coupling 
distribution. 

PACS: 75.10.Jm, 75.30.Ds, 75.50.Lk 
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I. INTRODUCTION 



The properties of quasi-particle excitations in low-dimensional disordered systems has 
been the subject of recent intensive investigations due to the possible applications of ran- 
dom structures and super-lattices in new devices. In many particular cases of interest, one- 
dimensional Hamiltonians can incorporate the main aspects of the disorder in what concerns 
the spatial distribution of quasi-particle excitations, such as optical, acoustic, electronic and 
spin-waves. Scaling theory establishes that in one-dimensional systems with uncorrelated 
disorder, the characteristic length of excitations is finite for any amount of disorder as a con- 
sequence of exponentially localized states [|TJ. Correlations in the disorder distribution can 
stabilize extended states. Resonant extended states emerge, for example, in one-dimensional 
electronic systems with randomly distributed impurity segments [0-^] and interpenatrating 
Anderson chains || . Recently, it has been demonstrated that long-range correlations in the 
disorder distribution can be responsible for the emergence of a phase of extended states 
within the band of allowed energies, with mobility edges separating localized and extended 
states ||[7|]. A recent optical experiment has demonstrated this phenomenon || which has 
been proposed for use in the development of window filters in electronic, acoustic, or photonic 
non-periodic structures 0. 

Spin-wave excitations in ferromagnetic chains with randomly distributed exchange cou- 
plings J have similar features as those of electronic excitations in chains with a particular 
distribution of pair-correlated off-diagonal disorder [|TU|-|il}]. The localization length of the 
low-energy one-magnon excitations diverges as £,(E) cx E~^, where <f) depends on the partic- 



ular form of the disorder distribution in the vicinity of J = pT| , |12[| . A diverging length of 
low-energy excitations is also observed in disordered harmonic chains [1^,16]. In both cases, 
this behavior is associated with the presence of an uniformly ordered ground state. 

Long-range correlations in electronic systems with off-diagonal disorder can stabilize ex- 



tended states in a more effective way than in systems with diagonal disorder fll7| . As the 
magnon equations of motion in ferromagnetic spin chains can be mapped onto those of 
electronic chains with off-diagonal disorder, one would expect the one-magnon excitations 
in disordered ferromagnetic chains to be also quite sensitive to the presence of long-range 
correlations in the random exchange distribution. In this paper we study, via direct diago- 
nalization and finite size scaling, the nature of the one-magnon eigenstates in a disordered 
chain with long-range correlated disorder with power spectrum decaying as S(k) cx k~ a . 
Our results indicate that a phase of extended magnons emerges at the bottom of the band 
for a > 1. We characterize the scaling aspects of the localization length and relate it with 
those of the disorder distribution. 



II. MODEL, FORMALISM AND NUMERICAL PROCEDURE 

The Hamiltonian describing a spin 1/2 quantum Heisenberg ferromagnetic chain with 
random exchange couplings has the form 

L 

h = — j n s„ • s n+ i , (i) 

n=l 
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the ground state of the chain being the ordered state with all spins aligned. The one-magnon 
excitations have the general form = J2n b n \n), where \n) represents a state with the spin 
on the n-th site flipped with respect to the ground state orientation. The coefficients b n can 
be shown to satisfy the difference equation 

{Jn-l + Jn)b n — Jn-lbn-l ~ J n b n+ i = Eb n (2) 

where E is the excitation energy. For an uncorrelated distribution of the exchange couplings, 
the one-magnon excitations are exponentially localized in the thermodynamic limite for any 
E > 0. The uniform E = mode is quite particular since this mode is not sensitive to the 
presence of disorder. This feature is responsible for the emergence of effectively extended 
low-energy excitations in finite chains |TT],|T2j. The typical localization length £ diverges as 
E —>■ with a power-law whose exponent depends on the specific form of the distribution 
P(J) at the vicinity of J = 0. For distributions with a finite average (1/J), the inverse 
localization length vanishes linearly with E. However, for distributions with a diverging 
(1/J), the excitations with low energy become more localized and only a slower divergence 
£ oc E~^ takes place with an exponent <p < 1 which depends of the specific behavior of P(J) 
at small values of J. This transition at E = 0, which exists even for uncorrelated disorder, 
is due to the fact that the ground state is not affected by the disorder. In contrast, this 
feature is not present in the usual electronic problem with uncorrelated disorder. 

In this work we will investigate the role played by long-range correlations in the exchange 
coupling distributions. In the analog electronic system, it has been recently demonstrated 
that long-range correlations can drastically modify the localized character of the excitations 
and even induce the emergence of a phase of delocalized states near the center of the band 

E3- 

In order to introduce long-range correlations in the disorder distribution, the exchange 
couplings J n will be considered to be in such a sequency to describe the trace of a fractional 
Brownian motion with a specified spectral density S(k) oc l/k a , where k is related to the 
wavelength A of the ondulations on the random exchange landscape by k = 1/A. For a = 0, 
one recovers the traditional disordered ferromagnetic chain model with 5-correlated disorder 
(J n Jn r ) = (J 2 )5 n ,n'- The exponent a describes the self-similar character of the random 
distribution and the persistent character of its increments. Following an approach based on 
the use of discrete Fourier transforms [p^ -p0|, a power-law spectral density is imposed by 
construction whenever the exchange couplings are given by 

L/2 
k=l 

where L is the number of sites and 4>k are L/2 random phases uniformly distributed in the 
interval [0, 27r] . The equation above is a general decomposition of a power-law correlated 
potential where randomness is present only in the relative phases of the Fourier components. 
We normalize the width of the disorder distribution to keep it finite and size independent 
with (J 2 ) — (J) 2 = 1. Such procedure implies that the couplings generated by Eq. 3 have to 
be rescaled by L~ a ^. In what follows, we shift the couplings to have (J n ) = 4 to enforce all 
generated couplings to be strictly positive. Those rare sequences not satisfying this bound 
where not considered on the subsequent analysis. The typical dependence of the coupling 
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constant landscape with the spectral density exponent a are similar to those reported in 
Ref. ||. The most relevant aspect is that the landscape becomes progressively less rough as 
a is increased, favoring therefore to derealization. 

Our numerical procedure consisted in obtain all eigen-energies and one-magnon eigen- 
states for finite chains of size L by direct diagonalization and to infer about the thermody- 
namic behavior through a finite size scaling analysis. We performed a configurational average 
over distinct disorder realizations such that, for each chain size, we computed 16 x 10 4 eigen- 
states. From each normalized eigen-state (J2 n b 2 n = 1), we computed the inverse participation 
ratio Y = J2 n ^t, an d used it to characterize the localized/delocalized nature of these one- 
magnon states. The typical localization length will be estimated by £ = l/Y. In general, 
the localization length is finite for exponentially localized states and diverges linearly with L 
for trully extended states. In the next section, we report our results for its scaling behavior 
for distinct long-range correlated coupling constants characterized by the power spectrum 
exponent a. 



III. RESULTS AND DISCUSSIONS 

The behavior of the inverse participation ratio Y and the typical localization length £ as 
a function of the one-magnon energy are shown in figure la and lb respectively as obtained 
from chains with L = 1600 sites and several power spectrum exponents a. For a = 0, 
which mimics an uncorrelated disorder distribution, we found that Y oc E at low-energy 
excitations, as expected. This trend remains in the range of a < 1. In this regime, we 
also found the localization length to be size independent over the entire band of allowed 
energy values (except for finite size scaling corrections near E — 0), even tough the local 
disorder scales down as L~ a ^. Our results indicate that the localization length of the 
excitations are not sensitive to such rescaling in this regime. In Fig. 2, we plot the normalized 
localization length ^{E)/L for a = 0.5 and several chain sizes to illustrate the comments 
above. £(E = 0)/L remains finite in the thermodynamic limite once the uniform mode 
remains delocalized. For any finite energy, the one-magnon excitations are exponentially 
localized with £ (E) / L — > as L increases. 

For a = 1, a new scaling behavior was identified where at low energies Y oc E^, with 
<p > 1, indicating that in this case the excitations have a stronger tendency to become 
delocalized. It is important to mention that a similar feature (0 > 1) was reported to be 
present in the analog electronic chain model with long-range correlated hopping amplitudes 
[L7|]. For any a > 1, the data suggest that a phase of delocalized low-energy magnon states 



emerges, which is seen quite clear in figure lb where £/L achieves a plateau. Therefore, long- 
range correlations can stabilize extended one-magnon states whenever the power spectrum 
exponent a > 1. 

To further characterize the scaling behavior of the localization length in the regime of 
strong correlations that allows for the existence of delocalized states, we employed a finite 
size scaling study of data from chains with L = 1600,800,400,200 and 100 sites. The 
behavior of £(E)/L is shown in figure 3a for the typical case of a = 2. For E < E c = 3.0(2) 
the data collapse characterizes a phase with delocalized states. E c depends on the correlation 
exponent a and on the disorder width. The one-magnon eigen-states remain exponentially 
localized for E > E c . However, the localization length in this regime becomes sensitive to 
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the size-dependent rescaling of the coupling constants. We found that the localized states 
for a = 2 have £ oc L 1 with 7 = 0.42 as illustrated by the collapse of data above E c in 
figure 3b. This size dependence is weaker than the expected for a similar rescaling of the 
local couplings in chains with uncorrelated disorder. The additional factor is due to the 
correlated nature of the disordered potential and its consequently enforced finite width. 
The new scaling exponent 7 depends on the correlation exponent a in a non-trivial way. 
In figure 4, we report our estimated values for 7. Our results suggest a non-linear relation 
7(a) oc (a — l) 75 , although we could not envisage a simple scaling argument to support 
this fact. In Figure 5, we show typical wave-functions representing the localized/delocalized 
nature of one-magnon excitations within each phase. 



IV. SUMMARY AND CONCLUSIONS 

In this work, we investigated the nature of one-magnon excitations in ferromagnetic 
quantum chains with long-range correlated disorder. We found that a phase of delocalized 
magnons emerges whenever the power spectrum exponent characterizing the strength of 
correlations is a > 1, in a close relationship with the behavior of the analog electronic 
model with off-diagonal disorder. Chains with distinct sizes and long-range correlations 
in the coupling exchange distribution need to have the local disorder scaled down as the 
size L increases in order to keep the disorder and band width finite. For a < 1, the 
localized spin excitations are insensitive to such local disorder rescaling. On the other side, 
for a > 1, the localized states acquire a size dependent localization length reflecting its 
sensitivity to the local disorder rescaling. However, such states should not be confused with 
critical multifractal states whose participation ratio also have a power-law size dependence 
2T|^24]. The present localized states have exponential tails and their anomalous scaling is 



linked to the scaling property of the underlying long-range correlated potential. It would be 
valuable to investigate the spin-wave dynamics in these new regimes and the possible relation 
between super- diffusion and correlation exponents. This information would be valuable for 
possible future spin-wave dynamics experiments on correlated ferromagnetic chains and/or 
non-periodic ferromagnetic super-lattices. 
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VI. FIGURE CAPTIONS 



Figure 1 - (a) The Inverse participation ratio as a function of the one-magnon energy E for 
a chain with L = 1600 spins. From top to bottom a = 0.0, 1.0, 1.5. The linear dependence 
Y oc E for low-energy excitations observed for a = is typical of uncorrelated potentials 
with finite < 1/J >. For a > 1 a phase of low-energy delocalized excitations emerges, 
(b) The scaled localization length £/L versus E for chains with L = 1600. From bottom to 
top a = 0.0, 1.0, 1.5, 2.0. The phase of delocalized states appears as a plateau for a > 1. 
Figure 2 The scaled localization length £/L versus E for a = 0.5. From top to bottom 
L = 100, 200, 400, 800, 1600. For excitations with finite energies, £(L)/L — > as L increases, 
characterizing localized states. The uniform mode at E — remains delocalized. 
Figure 3 - (a) The scaled localization length £/L versus for a = 2.0. From top to bottom 
L = 100,200,400,800, 1600. The phase of delocalized states appears as a size independent 
plateau. The critical energy separating localized and delocalized states is E c = 3.0 ± 0.2. 
(b) The scaled localization length £/L ' 42 versus E for a = 2.0 showing a data collapse for 
E > E c . From bottom to top L = 100,200,400,800,1600. The size dependence of the 
characteristic length of such exponentially localized states reflects their sensitivity to the 
rescaling of the local disorder. 

Figure 4 - Estimated values of the exponent 7 related to the anomalous scaling of the 
localization lenght above E c for several values of the correlation exponent a. The data 
suggests 7(a) oc (a — l) a75 . 

Figure 5 - Typical wave-functions representing delocalized and localized one-magnon states 
of a long-range correlated ferromagnetic chain with a = 1.75. (a) Delocalized state near 
E = 2.0; (b) Localized state near E = 8.0. 
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